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A non-minimal coupling between the dark matter and dark energy components may offer a way of 
solving the so-called coincidence problem. In this paper we propose a low -2 test for such hypothesis 
using measurements of the gas mass fraction /gas in relaxed and massive galaxy clusters. The test 
applies to any model whose dilution of dark matter is modified with respect to the standard a~^ 
scaling, as usual in interacting models, where a is the cosmological scale factor. We apply the test to 
current /gas data and perform Monte Carlo simulations to forecast the necessary improvements in 
number and accuracy of upcoming observations to detect a possible interaction in the cosmological 
dark sector. Our results show that improvements in the present relative error (Tgas/Zgas are more 
effective to achieve this goal than an increase in the size of the /gas sample. 

PACS numbers: 98.80.-k; 95.35.-|-d; 95.36.+x; 98.65.Cw, 98.80.Es 


I. INTRODUCTION 

One of the most important goals of current cosmolog¬ 
ical studies is to unveil the nature of the mechanism be¬ 
hind cosmic acceleration. In principle, the most promis¬ 
ing explanations involve either the existence of new fields 
in high energy physics or large-scale modifications of 
gravity theory. Following the former route, several possi¬ 
bilities have been discussed in the literature, among them 
a possible interaction between the dark matter and the 
smooth dark energy component, permitting transfer of 
energy between these two fields and violating adiabatic- 
ity (see [T] and references therein). 

If the dark energy component is the vacuum energy 
density p\, characterised by an equation-of-state param¬ 
eter w = p/p = —1, where p is the pressure, these sce¬ 
narios may constitute not only a phenomenological at¬ 
tempt at alleviating the so-called coincidence problem, 
since dark matter and dark energy are no longer assumed 
to be independent, but also the unsettled situation in the 
particle physics/cosmology interface, in which the cos¬ 
mological upper bound for pA differs from theoretical ex¬ 
pectations (p„ ^ lO^^GeV'^) by more than 100 orders 
of magnitude (see [5] for a discussion on the cosmolog¬ 
ical constant problems). From the physical viewpoint, 
it worth mentioning that no known symmetry in Nature 
prevents a non-minimal coupling between these compo¬ 
nents. Therefore, although no observational piece of ev¬ 
idence has so far been unambiguously presented for an 
interaction in the cosmological dark sector, a weak cou¬ 
pling still below current detection cannot be completely 
excluded. 

Given the theoretical implications of such a detec¬ 
tion for cosmology, it is important to investigate ways 
to probe this possibility from observations. Our goal in 
this paper is to propose and explore a low-z test for a 
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possible dark matter/dark energy interaction from mea¬ 
surements of the gas mass fraction /gas in relaxed and 
massive galaxy clusters. The test proposed is model- 
independent in the sense that in the low-z regime the /gas 
equations become independent of the interacting model, 
although the signature of a non-minimally coupled re¬ 
mains present. We apply the test to the current data 
and investigate the constraining power of future /gas mea¬ 
surements from Monte Garlo (MG) simulations. We also 
discuss the necessary conditions in terms of number and 
accuracy of upcoming data for a successful application of 
the test. In what follows, we outline the main assump¬ 
tions of our analysis and discuss the main results. 

II. INTERACTION IN THE DARK SECTOR 

In ou analysis, we assume a non-minimal coupling be¬ 
tween the cold dark matter (CDM) and a dark energy 
component (A) with w = —1 (baryons and radiation are 
assumed to be separately conserved, implying the usual 
scaling pb oc a~^ and p-^ oc a~^). The evolution of the 
mean dark matter density pm and of pA are obtained 
from the balance equations 

Pa = -Q, (1) 

Pm + ^Hpm = —Q , (2) 

where a dot denotes derivative with respect to time, 
H = ^ is the Hubble parameter and a is the cosmo¬ 
logical scale factor, written in terms of the redshift pa¬ 
rameter ds a = 1/{1 -\- z). In the absence of a natural 
guidance from fundamental physics on the interaction 
function Q, which controls the rate of energy transfer, 
several phenomenological forms have been proposed [T]. 
An interesting approach was discussed in Ref. [3] in which 
the time variation of A is deduced from its effect on the 
CDM evolution. The qualitative argument is the follow¬ 
ing: since A is decaying into CDM particles, CDM will 
necessarily dilute more slowly compared to its standard 
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evolution, pm oc a Thus, if the deviation from the 
standard evolution is characterised by a constant e, i.e., 

Pm oc (p{a)a~^ where ^p{a) = , (3) 

Eqs.Q and ([^ yield 

PA = A — -a ^ + const. , (4) 

O 6 

which is equivalent to the dark energy evolution found 
if one assumes Q = epmH (from now on, the subscript 
0 denotes present-day quantities). Note that a positive 
(negative) value of e leads to a slower (faster) dilution of 
the matter component, which implies in an energy trans¬ 
fer from A (CDM) to CDM (A) (for thermodynamic con¬ 
straints on the sign of e, we refer the reader to mm)- 
A more general parameterisation for the interaction pa¬ 
rameter was discussed in Ref. [5], in which this quantity 
is assumed to be an arbitrary function of the scale fac¬ 
tor, i.e., e(a). In this case, regardless of the interaction 
function Q a given theory may provide or, equivalently, 
of its effect on the CDM evolution, pm can always be 
parameterised in terms of e(a) = log^^ (p(a). Note also 
that this approach can be applied even if the dark energy 
component has an arbitrary equation-of-state parameter 
(w(z) ^ -1). 


III. A LOW-z TEST FOR INTERACTING DARK 
ENERGY 


have been adopted in the literature (see, e.g., [Sl ITTHT^ l. 
However, for the test discussed here the exact form or 
value of C/ is not important, so that we marginalise over 
it in all the analyses presented in the next sections. 

Regardless of the cosmological models involved in 
the analysis, the distance ratio above equals unity 
in the cosmographic limit (z << 1). Moreover, 

if e = 0, the baryon-to-total-mass density ratio 

^b{z)/[^b{z) -I- ^lm{z)\ becomes constant, as the baryons 
and CDM density parameters have the same 2 ; depen¬ 
dence. On the other hand, under the assumptions dis¬ 
cussed earlier, if there is a non-minimal coupling in 
the cosmological dark sector the above equation can be 
rewritten as 


fgas{z << 1 ) 





( 6 ) 


where the extra term {l + z)~'^^^'> makes fgas{z « 1) no 
longer constant, providing a direct test to the idea of an 
interacting dark energy field. In reality, this test can be 
applied to any scenario whose dilution of dark matter is 
modified with respect to the usual a~^ scaling, regardless 
whether it is due to an interacting or any other physical 
process. Here, however, for simplicity, we will restrict our 
analysis to the case in which a non-minimal coupling in 
the dark sector is parameterised by a constant interacting 
parameter e, as given in Eq. (§. 


Observations of the gas mass fraction in relaxed and 
massive galaxy clusters have been widely used as a cos¬ 
mological test (see, e.g., [MS] and references therein). 
The key assumption of this test is that the baryonic mass 
fraction in galaxy clusters provides a fair sample of the 
distribution of baryons at large scale M- The mass frac¬ 
tion of hot gas measured within a characteristic radius of 
a halo at redshift z can be written as 

3 
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IV. IMPLEMENTATION OF THE TEST 

A. Data 

In our analysis, we consider the Chandra fgas data 
analysed in m and shown in Fig. la. Originally, this 
sample comprises 42 regular, relatively relaxed and mas¬ 
sive galaxy clusters distributed over a wide range of red- 
shift (0.063 < z < 1.063). 


B. Gurrent constraints 


fgas{z) = G 


H{,(z) 


^h(z) -f Glm.iz') 


where Of,(z) and flm(z) are, respectively, the baryon and 
CDM density parameters as function of z and the an¬ 
gular diameter distance (cIa) ratio accounts for devia¬ 
tions in the geometry of the universe between the hdu- 
cial model adopted in the observations and the model 
to be testecQ The term G accounts for star formation 
and other baryon effects within the characteristic radius 
under the assumption that fgas should be approximately 
constant with redshift (for a discussion of the systemat- 
ics errors included on G, see also m)- Many forms of Q 


^ For a general derivation of fgas{z) taking into account possible 
violations of the so-called distance duality relation, see m- 


In order to implement the low-z test of Eq. (§, 
we first study the redshift interval in which the cosmo¬ 
graphic condition above is satisfied. The fiducial model 
used in the observations was a flat ACDM scenario with 
m- Thus, considering Gib = 0.0489, given by 
current CMB experiments m, we show in figure lb the 
quantity <5 = [(df X 100 %, where 

eCDM stands for the interacting model discussed in Sec. 
H with the values of e and Glm,o fixed at the best-fit values 
given by Eq. We consider fgas data up to z = 0.4, 
which corresponds to a dispersion of 0.6% for the angu¬ 
lar diameter ratio of Eq. (§. Such value is negligible 
when compared with current observational uncertainties 
fgas/fgas — 11.4%) and restricts our working sample 
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FIG. 1. Left) Gas mass fraction measurements as a function of redshift for 42 galaxy clusters studied in Ref. [TT]. Solid circles 
corresponds to the reduced sample (27 data points) obtained after applying the cosmographic condition discussed in the text. 
Middle) Dispersion of the angular diameter distance as function of the redshift. The maximum limit allowed is <5 = 0.6%, which 
reduces our sample to clusters up to z < 0.4. Right) Evolution of the figure of merit defined in the text as a function of the 
number of data points N. The current value, obtained from the 27 data points shown in the panel of the left, is FoMoSs — 425 


to 27 clusters (Fig. la)|^ 

Considering this sub-sample, we perform a analysis: 


27 




where is the theoretical expression for fgas at low- 
z given by Eq. (j^ and is the error in the 

measurements. From this analysis, we find 


= 0.28 ±0.01 and e = 0.035+^;235 (Icr) , (8) 


with xt = 1-17 stands for the number of degrees of 
freedom). Clearly, the current data are unable to unam¬ 
biguously detect or rule out a non-minimal coupling in 
the cosmological dark sector. In what follows, we per¬ 
form MC simulations to discuss the necessary conditions 
(i.e., number of fgas data points and improvements on the 
relative error measurement crfgas) to detect such a 
coupling at a given confidence level. 


V. SIMULATIONS AND RESULTS 

We expect that the larger the number of fgas data 
points or the smaller their error the smaller the 


^ Note that this large z interval is possible only because the fidu¬ 
cial model used in the fgas observations and the best-fit scenario 
found below [Eq. predict very similar dj^{z) evolution. How¬ 
ever, the most general way to define the test of Eq. is through 
the cosmographic limit discussed above. 


uncertainty on the interacting parameter should be. 
Thus, in order to find which are the necessary improve¬ 
ments in future observations to check the hypothesis of a 
non-minimal coupling in the cosmological dark sector, we 
generate synthetic samples of fgas measurements assum¬ 
ing the distribution and redshift interval of the current 
observations. 

The steps of the simulations are the following: we first 
define the number of points {N) to be simulated. Then, 
we simulate a sample of fgas{z) measurements, with N 
points in the same redshift range of the current data. 
For each point simulated we assume a normal distribu¬ 
tion a), with a mean value and error defined as the 
fiducial values. The relative errors for fgas are obtained 
assuming that the relative error provided by observations 
can be fitted by a polynomial to find its redshift depen¬ 
dence. We performed 10^ realisations for each one of the 
analyses described below. 

We first study the evolution of the Figure-of-Merit, 
defined as FoM = as function of the number 

of data points in the range [50 - 250]. Using the simulated 
data we perform x^ analyses and find the best-fit values 
for the pair (e, We repeat this process 10^ times 

and calculate the mean and standard deviation of these 
parameters and the corresponding FoM, for a fixed value 
of N. The fiducial model used in the simulations has the 
parameters values given in Eq. Q and the results are 
shown in Eig. Ic. From this analysis, we find that for a 
relatively small number of measurements when compared 
to future X-ray survey^ e.g., N = 200, the FoM value 


® See, e.g., http://www.mpe.mpg.de/eROSITA 
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FIG. 2. Left The interacting parameter e as a function of the number of data points N. Assuming the values given in Eq. Q 
as fiducial value, the number of points for which a,, = e is Ncrit = 440. Right Percentage of the relative error as a function of 
the critical number of data points. For crf^^^/fgas = 100%, the results of the previous panel are recovered. 


1.5 X 10^ corresponds to an improvement of order of 3 
with respect to the current data (FoMobs — 425). This 
clearly shows that very tight bounds on the interacting 
parameter can be obtained from a larger sample of galaxy 
clusters X-ray observations. 

In what follows, we will direct our analysis to the neces¬ 
sary improvements in future observations, i.e., the num¬ 
ber of galaxy clusters N and the relative error of the gas 
mass fraction <Jf^^Jfgas, to observationally distinguish 
between values of e of the order of current estimates and 
the standard value e = 0. Besides assuming the best-fit 
values given in Eq. Q as fiducial values, our simulations 
always assume the distribution and relative error of the 
current fgas sample shown in Fig. la. 

Under such conditions, the simulation process provides 
the number of data points N^rit necessary to verify the 
hypothesis of a non-minimal coupling in the cosmological 
dark sector. In order to find Ncrit, we adopt the criterion 
< 1, where n corresponds to the confidence level 
(number of standard deviations) required. The results 
are shown in the Fig. 2a. As expected, the error bar for 
the interaction parameter decreases for higher values of 
N. For n = 1 (Icr), we find Ncrit = 440, which clearly 
illustrates the viability of implementation of the test pro¬ 
cedure in the near future since a value of >> Ncrit 
is expected in some planned X-ray surveys. It is worth 
noticing that such result is somehow conservative since no 
improvement on the upcoming fgaa data relative to the 
current ones was considered to derive the above value of 

Ncrit• 

In order to study the influence of the uncertainties of 


fgas measurements on the estimates of Ncrit, we also 
vary the relative error fgas with A, considering 

the interval between 50% — 100% of the current obser¬ 
vational value. This is shown in Fig. 2b, where 100% 
means the present-day value. Clearly, by improving the 
quality of the data or, equivalently, decreasing their rel¬ 
ative error, Ncrit becomes considerably smaller. As an 
example, decreasing the relative error in 50%, we find 
Ncrit = 145, which corresponds to 1/3 of the previ¬ 
ous value obtained when the relative error of the cur¬ 
rent fgas data fgas = 100%) is considered. For 

Ncrit = 27, which corresponds to the number of fgas 
measurements currently available, the improvement on 
^fgas/ fgas needed to distinguish coupled from uncoupled 
dark energy is of the order of 25% relative to the present 
observational value. These results indicate that improve¬ 
ments in the relative error of fgas observations are more 
effective to achieve this goal than an increase in the size 
of the sample. 


VI. CONCLUDING REMARKS 

Testing a possible non-minimal coupling between the 
dark matter and dark energy fields constitutes an im¬ 
portant task for cosmology and fundamental physics. 
In this paper, we have proposed and applied a new, 
model-independent test for this hypothesis using low-x 
measurements of the gas mass fraction in galaxy clus¬ 
ters. As well known, the transfer of energy between the 
dark components necessarily modifies the dilution of dark 
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matter with respect to the usual a~^ scaling whereas 
the baryonic component remains separately conserved 
{pb oc a~^). Therefore, the test here proposed can in 
principle be applied to any model of this class or even to 
cosmologies in which the attenuated dilution of the CDM 
component results from another physical process. This 
is possible because in the low-z regime the /gas equa¬ 
tions become cosmological model-independent, although 
the signature of a non-standard CDM evolution remains 
present. 

Since the current observational data cannot unambigu¬ 
ously detect any interaction between the dark compo¬ 
nents, we have used the present bounds given by Eq. Q 
as fiducial values and performed MC simulations to study 
the necessary requirements for a e 7 ^ 0 detection. We 
have then generated thousands of /gas samples assuming 
the present observational error distribution, as described 
in Sec. V. After discussing the gain in the FoM with 
the number of data points N in the synthetic sample, we 
have investigated the number and accuracy of future ob¬ 
servations needed for a successful application of the test. 
For the simplest case in which e = const., adopted in our 
analysis, and considering the relative error Of^^^/fgas of 
current data, we have found that the minimum number of 
data points necessary to distinguish between an interact¬ 


ing and a non-interacting dark energy component should 
be ^ 15 times larger than the current observational sam¬ 
ple {Ncrit = 440 and A/hs = 27). This number, how¬ 
ever, is much smaller than the number of galaxy clusters 
expected to be observed by the extended Roentgen Sur¬ 
vey with an Imaging Telescope Array (eROSITA), which 
clearly shows the viability of implementation of this fgas 
test in the near future. We have also calculated how 
Ncrit changes as the relative error of fgas measurements 
improves. As shown in Fig. 2b, decreasing fgas by 

a factor of 2 leads to Writ = 145, which is three times 
smaller than the value found assuming the relative error 
of the current fgas data. These results clearly show that 
a combination of quantity and precision of future fgas 
measurements can become an important tool not only 
to probe a possible interaction in the cosmological dark 
sector but also to explore its consequences. 
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